Anisotropic Ginzburg-Landau and Lawrence-Doniach Models for 

Layered Ultracold Fermi Gases 

Mauro lazzi/ Stefano Fantoni,^ and Andrea Trombettoni^' ^ 

^SISSA, Via Bonomea 265, 1-34136 Trieste, Italy 
^ANVUR, P.za Kennedy 20, I-OOlU Roma, Raly 
^CNR-IOM DEMOCRITOS Simulation Center, 
Via Bonomea 265, 1-34136 Trieste, Raly 

We study the anisotropic Ginzburg-Landau and Lawrence-Doniach models describing 
a layered superfluid ultracold Fermi gas in optical lattices. We derive the coefficients 
of the anisotropic Ginzburg-Landau and the mass tensor as a function of anisotropy, 
filling and interaction, showing that near the unitary limit the effective anisotropy 
of the masses is significantly reduced. The anisotropy parameter is shown to vary 
in realistic setups in a wide range of values. We also derive the Lawrence-Doniach 
model - often used to describe the 2D-3D dimensional crossover in layered super- 
conductors - for a layered ultracold Fermi gas, obtaining a relation between the 
interlayer Josephson couplings and the Ginzburg-Landau masses. Comparing to the 
Ginzburg-Landau description, we find that the region of validity of the Lawrence- 
Doniach model is near the unitary limit. 
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The presence and nature of the superfluid and superconducting transition in systems 
of different dimensionahty is a widely studied and interesting topic. It is well known that 
the transition in two and three dimensions is radically different, as the first is governed by 
a Kosterlitz-Thouless mechanism of topological origin (the unbinding of vortices) and the 
second is the paradigmatic example of the symmetry breaking well described, at least qual- 
itatively, by mean-field theories. This is the main reason of the wide interest in the study of 
materials with properties in-between the two, i.e. layered superfluids and superconductors. 
Their phenomenological properties are usually well described by anisotropic effective mod- 
els: a successful example is given by the Lawrence-Doniach (LD) model |Tj (in which the 
discrete layers are weakly coupled through Josephson terms) describing the electrodynamics 
of layered superconductors and copper oxides [2] . Another related way to introduce the ef- 
fect of the layering is to consider a Ginzburg-Landau (GL) theory having anisotropic masses 
(eventually phenomenologically determined) [2]. A well known feature of these phenomel- 
ogical models is the fact that they work well even if the underlying microscopic mechanisms 
are not known (as it occurs for high-Tc cuprates). 

In this paper we study and compare the anisotropic GL and LD models describing a 
layered superfluid ultracold Fermi gas in optical lattices: the coefficient of these phenomeno- 
logical models can be microscopically derived and the parameters, like the strength of the 
optical lattices, can be tuned with high accuracy in realistic experimental setups. We deter- 
mine the mass tensor in the GL theory as a function of anisotropy, filling and interaction, 
showing that near the unitary limit the effective anisotropy of the masses is significantly 
reduced. By tuning the strength of the optical lattices, the dimensionless anisotropy param- 
eter 7 in realistic setups can vary in a wide range of values: 7 ~ 1 — 200 (for comparison 
with layered superconductors, we remind that 7 ~ 7 for YBCO and 7 > 150 for BSCCO). 
We also derive the LD model, discussing a relation between the interlayer Josephson energy 
couplings in the LD model and the masses in the GL model. Using this result, the LD model 
with only couplings between adjacent planes is shown to be equivalent to the anisotropic 
GL theory in the BEG limit, while in the BGS side (where the size of the pairs is larger than 
the lattice spacing) next-nearest-neighbour couplings have to be taken into account. 

Layered superconductors and superfluids: The first example of a layered superconductor 
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was found in a naturally occurring compound of transition-metal dichalcogenides layers 
(such as TaS2) intercalated with organic molecules |3], where, for temperatures of a few 
Kelvin, the metallic planes become superconductors. The discovery was followed later by 
the creation of artificial samples with alternating layers of different metals, having a different 
transition temperature 01 . The nature of the transition was found to be three dimensional 
[5], and the system was shown to be very well described by an anisotropic GL free energy 
at and above the transition, where the coherence length diverges as ^ oc (T — Tc)^^/^. At 
lower temperatures where the length ^ becomes comparable with the plane separation, the 
presence of a layer structure becomes important. To describe this Lawrence and Doniach 
proposed a tight-binding approximation [1] able to describe the crossover between the 2D 
behaviours of essentially decoupled planes to the 3D. 

Another very important class of superconductors having a layered structure is the one 
of cuprates, discovered in 1986 by Bednorz and Miiller [6]. Their characteristic structure 
presents a layer of one or more conducting copper oxide planes alternated to an insulat- 
ing lattice with a perovskite-like structure. The layered structure makes the behaviour of 
cuprates strongly anisotropic and the dimensional crossover temperature is in general very 
near the critical temperature for the resistive transition [2]: despite the lack of a definitive 
consensus on the mechanism responsible for the superconductivity in the cuprates, it is by 
now established that the LD model successfully describe the transition and the electromag- 
netic response of cuprates |2l El E] . 

The study of the dimensional crossover in cuprates is not in general an easy task: on one 
hand, the smallest barrier attainable is that of a single plane of an insulating material, on 
the other hand inserting several planes is a delicate process as the separation can become 
nonuniform [8]. 

From this point of view, ultracold atoms (and in particular ultracold fermions) in optical 
lattices provide a favourable setting: the tunnelling rate t can be controlled by adjusting 
the optical lattice parameters [9l [10] and the on-site energy coefficient U can be controlled 
through Feshbach resonances |TT1 [12]. The ratio U/t can be then controlled with high 
precision, as well as the geometry and the anisotropy of the system (e.g., optical lattices with 
different strengths can be superimposed in different directions). For ultracold bosons, optical 
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lattices has been used to induce a Mott-insulator/superfluid quantum phase transition [13] 
and to study Josephson dynamics in many- and two- wells settings [El [15]. The study of 
ultracold fermions in optical lattices also stimulated a huge amount of experimental work 
[T6H2T] and they provide a field of research in which even more far-reaching experiments 
are expected to take place. The low energy properties of ultracold fermions in deep optical 
lattices are described by Hubbard-like models [22]: metallic and insulating phases in 3D 
lattices have been studied [HI [20]. in presence of attractive interaction {U < 0) one has 
then a physical realization of the attractive Hubbard model [231427] . 

Superimposing a ID optical lattice to an ultracold Fermi gas gives raise to a stack of 
two-dimensional layers: very recently, the evolution of fermion pairing from 3D to 2D was 
experimentally studied, showing the opening of a gap in the radiofrequency spectrum when 
the dimensionality is reduced ^21j. The confinement of Fermi mixtures in 2D was also ex- 
perimentally investigated [281433] . showing properties different from the 3D case [31]. The 
observation of a pairing pseudogap in a 2D Fermi gas, using momentum-resolved photoemis- 
sion spectroscopy, has been reported [33]. An important advantage of the use of an optical 
lattice is that when it is applied (say along the z direction) then the interlayer coupling 
between 2D pancakes along z can be easily modified by varying the strength of the optical 
lattice. Furthermore, if two additional optical lattices are added in the xy plane, then the 
tunnelling t± can be tuned separately from t\\ and the layering can be tuned by varying the 
ratio t\\/t±. 

Effective models: Near the critical temperature, the superfiuid wavefunction ^ is small, so 
that the free energy difference 5F between the normal and superfiuid state can be expanded 
in powers of [35]: for a layered superconductor /superfiuid 5F reads then 

= I rff |^(|9.vl>|2 + |5^vl>|2) + ^|9,vl>|2 + a{Tm' + /3|v]>|4| . (1) 

The phenomenological coefficients a and /3 appearing in ([T]) can be derived from experiments 
or from the microscopic Hamiltonian; the masses are those of the Cooper pairs. The 
GL equations describing the equilibrium state can be found by minimizing ([T]) with respect 
to \E'. For a > there is only the minimum \l/ = corresponding to the normal state, 
while for a < there are uniform solutions = |\E'o|e*^, hence the critical temperature is 
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given by a{Tc) = 0. In the LD tight-binding approach one has a 2D superfluid wavefunction 
for in the nth plane along z and the term with transverse derivative becomes 
""'^i l^n(^) ~ ^n+i('")P with the nearest-neighbour tunnelling coefficient Ji (xl/d? {d is 
the interlayer separation, i.e. the lattice spacing of the periodic potential along z). 

In layered superconductors, the effects and the penetration of a magnetic field have been 
thoroughly studied [2[[7l|36]: such effects can be incorporated via the minimal substitution 
di di — i^Ai (where A is the vector potential) and adding to the free energy the mag- 
netization energy Fh = {B — H)^/87i. The penetration of the magnetic field in anisotropic 
type-II superconductors and the response of a layered superconductor to an electromagnetic 
field have been an established field of research for decades [21 [3 |36]. In the following we 
show that the layered GL theory ([T| holds for a layered ultracold Fermi superfluid gas: then, 
to reproduce the phenomenology of the electrodynamics of layered superconductors in the 
context of ultracold atoms one can use a layered Fermi mixture in an fictitious magnetic 
field realized in a synthetic way using spatially dependent optical couplings between internal 
states of the atoms f37J or by rotating the trapping potential [38] . 

The Hamiltonian describing a two-component Fermi gas in a layered optical lattice reads 
H — fiN, with the Hamiltonian H given by 

H = -Y. + h.c.)-uY, ^lA^ikv (2) 

{ij)<7 i 

where N = J^ia^lo-^ia- the number operator, /i is the chemical potential, i,j are site 
indices, a = t) i is the spin index and is the fermionic operator (the sum is on distinct 
pairs of nearest-neighbours). The interaction is assumed attractive, U > 0, and the total 
number of fermions per site is denoted by n. The physical system modelled by ^ considers 
two optical lattices in the x, y directions (determining the tunnelling rate between the 
sites in the xy planes) and a different lattice in the z direction, with interlayer coupling t^. 
Notice that in our subsequent determination of the coefficients of the anisotropic GL and 
LD models the crucial ingredient is to have the optical lattice along z and the mechanism 
of layering in the form of tunnelling between adjacent layers: if the optical lattices in the 
x,y directions are absent (or small), then the anisotropic GL ([T| is again obtained (with 
my = 2m, where m is the mass of the fermionic atoms). Similarly, in writing ^ we assumed 
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FIG. 1. Critical temperature ksTc as function of U (both scaled in units of the total bandwidth 
D) at various particle densities - n = 0.9 (solid lines), n = 0.5 (dashed lines), n = 0.1 (dotted lines) 
- and anisotropics tx/ty = 1, 0.3, 0.1 (in each set from top to bottom). 

that in each minimum of the lattice there is a total number of fermions n < 2 (i.e., only the 
lowest band of the periodic potential is occupied): otherwise, we would have a multi-band 
Hubbard model and again ([T]) would be retrieved. 

The single-particle energies of Hamiltonian for U = are = — 2t|| (cos k^d + 
cos kyd) — 2t± cos kzd with total bandwidth D = Sty +4t_L. The isotropic 3D case corresponds 
to t|| = t_L, while the isotropic 2D case to t_L = (isolated pancakes). Estimates of the value 
of U for which the BCS-BEC crossover occurs gives U ~ 0.7D in 3D and U ~ 0.2D in 2D 
[23H27]. 

To derive the coefficients of the anisotropic GL we proceed in a path-integral setting 
m] : performing the Hubbard-Stratonovich transformation and integrating out the fermionic 
field, a bosonic effective action is obtained (the critical temperature is the one at which 
the order parameter is vanishing). Details are in Appendix A: the coefficients of the GL free 
energy ([T]) are found to be 
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X, z denotes the direction and the mass parameters Ma are related to the masses 



vfiy = my, niz 



m±). It is e-j. 



in the GL energy Q by Ma = maW^ (with m, 
X = tanh(/3ceg/2) and Y = cosh-2(/3ceg/2). 

It is also possible to derive the LD tunnelling coefficients: in the LD model, the xy layers 
are weakly coupled along the z direction through Josephson terms. Usually, only Joseph- 
son couplings between adjacent planes are considered. Here we consider the more general 
possibility of longer range tunnelling: this amounts to consider a Josephson contribution 
~ "^".m I '^n — ipm P to the energy, where n,m are denoting the different planes (i.e., 
the minima along z). We find 



J n 



27r 



(lq^(,-iin-m)q-.d^^-^ = 0; 0, 0, q,) 



(5) 



where the expression for x(g) = x(^i^; (f) is reported in Appendix A, Eq.(AlO). In presence 
of lattices in the x,y directions, one can similarly define Josephson couplings jJl.m. It is 
possible to determine a useful relation between the coefficients Jn,m and the mass m_L of the 
anisotropic GL model: one finds 



2m±d'^ 



(6) 



p=i 



where we used the notation Jp = Jn,n+p (we are assuming translational invariance along z): 
then Ji = Jn,n+i- 
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FIG. 2. Anisotropy parameter 7^ = m±/m\\ vs. U/D for different anisotropies: from top to bottom 
t±/ti\ = 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 0.7, 0.8, 0.9, 1 (with n = 0.5). 
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Results: The quantities X, Y entering Eqs. (|3])-(|4]) depend on (3c = l/fc^Tc: therefore 
in order to evaluate the masses my, m_L one has to determine Tc and as a function of 
the interaction U, the filhng n and the anisotropy t^\/t±. In Appendix A we write the two 
equations from which we can determine Tc and fi. The resuhs for the critical temperature 
Tc for various values of filling n and of the ratio are plotted in Figjlj we plot both Tc 
and U in units of D, showing that in these units the effect of the anisotropy is relatively 
small (the critical temperature in the interval 0.5 < t_L/t|| < 1 is very similar to the Tc of 
the isotropic case). The maximum value of the critical temperature occurs practically at 
the same value of the interaction when the total bandwidth D is used as energy unit. 

Using the critical temperature Tc, we can compute the mass anisotropy parameter 7 
usually defined [2] as 

2 _ (7) 



7 



mil 



The anisotropy parameter 7 plays a crucial role in the electrodynamics of layered supercon- 
ductors since a formula giving the value of a quantity for a general anisotropy parameter 7 
given its value for 7 = 1 holds [42j . Values of 7 for several superconducting compounds has 
been determined from experimental data j2]. In layered Fermi superfluids, the value of 7 
affects the critical value of a synthetic magnetic field (this is the atomic counterpart of the ef- 
fect of a magnetic field on a layered superconductor [2j): supposing to have a general rotation 
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FIG. 3. Plot of the quantity 2'm±cPJi/h^ vs. U/D at n = 0.5 (with d lattice spacing) for the same 
values of t±/t\\ of Figj2| 
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fi, the critical values of the rotation for the vortex melting obey the relation f2c2||/^c2± = 7- 
The mass ratio could be experimentally determined from momentum-resolved photoemis- 
sion spectroscopy [33] or from the measurement of the frequency of the collective excitations 

m- 

The results for the anisotropy parameter 7 of a layered ultracold Fermi gas are shown in 
Figj2] 7 significantly decreases close to the unitary limit due to the effect of the increase of 
Tc in the crossover region. This result is consistent with the qualitative expectation that the 
effect of external anisotropics is softened in the unitary limit. Results for 7 vs. filling are 
reported in Fig|5| the filling modifies 7 on the BCS side, while the anisotropy on the BEC 
side remains virtually identical. To understand if and how much the collective behaviour 
enhances or not the (external) anisotropy artificially tuned by imposing a layering through a 
value of smaller than one, we studied an effective anisotropy defined as 7^^^ = 7^-t±/t||. 
The results are drawn in Fig|6] it is seen that in the BCS side the enhancement is more 
pronounced with respect to the BEC side. Again, the effective enhancement is reduced in 
the unitary limit. We also checked that in the BEC limit the quantity 7^ ■ (t_L/t||)^ is found 
to approach 1 for different fillings, see FigjT} 
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FIG. 4. Plot of the validity region of the LD model (with n = 0.5), obtained by fixing an arbitrary 
value Jf = 0.8 ^^'^2 at which we consider the Josephson tunnelling to be a good description of the 
pair motion along axis a = x,y,z (here Jf = Jf = j| and Jf = Jj*"). The shaded region is where 
the pairs move with discrete jumps in the z direction, but not in the x and y directions. 
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Finally, we numerically evaluated the LD coefficients (|5]) and compared them to the GL 
masses in Eq. ([s]). We found that in the BEG side the nearest- neighbour Josephson coupling 
Ji is enough to give the correct value m_L, but not on the BGS side of the crossover, where the 
size of the pairs is significantly larger than d. Then, the anisotropic GL model is equivalent 
to the LD equations with the Josephson couplings Ji only in the BEG regime (for arbitrary 
filling): using only nearest-neighbour Josephson couplings for the free energy instead of the 
continuous kinetic energy term in the direction z is approximately correct when 2mj^d?Ji/f? 
is near 1, as one can see from Eq. ([6]). The results in Figjs] show a non- monotonous 
behaviour around the unitary limit. The validity region of the LD model with only Josephson 
couplings Ji along z is seen to increase when the layering becomes stronger (i.e., t±/t\\ 
becomes smaller). In presence of lattices in the x,y directions, the use of the LD model 
is especially appropriate when the Josephson couplings can be used in the perpendicular 
directions (across the planes) but not in the parallel direction (between the planes). To give 
a qualitative idea of the interval where the LD model is valid one can choose an arbitrary 
threshold (we choose 0.8) and plot the region where 2m\\<fjl/h^ < 0.8 < 2m^<fJi/h^. As 
can be seen in Fig|4| the validity region of the LD model is narrower for small layering (i.e., 
t±/t\\ close to 1) and it is centered around the interval in which the critical temperature 
reaches its maximum value. This result is pictorially summarized in Figjs] 

Conclusions: In this paper we studied the superfluid phase of an ultracold layered Fermi 
gas in the presence of an optical lattice near the critical temperature. This system is de- 
scribed by an anisotropic Ginzburg-Landau theory (that is also used for other systems such 
as layered superconductors). We derived the coefficients of the Ginzburg-Landau equation: 
the anisotropy was shown to be enhanced by the superfluid, and the anisotropy parameter 
7 decreases near the crossover region. The dependence on filling was found to be important 
on the BGS side of the crossover and negligible when the gas is in the molecular condensate 
phase. These results link the underlying microscopic structure to the macroscopic proper- 
ties of the superfluid: it is then possible to tune the experimentally controllable parameters 
such as interaction strength, filling and lattice strength in order to emulate the phenomeno- 
logical electrodynamics of layered superconductors. We also derived the Lawrence-Doniach 
model for the layered ultracold Fermi gas and we discussed a relation between the interlayer 
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Josephson couplings in the Lawrence-Doniach model and the masses in the Ginzburg-Landau 
energy. Using this result, the Lawrence-Doniach model with only couplings between adjacent 
planes is shown to be valid in the crossover region between the BCS and the BEG limits. 

Acknowledgments: Discussions with K. Schmidt and G. Roati are gratefully acknowl- 
edged. 



Appendix A: Coefficients of the anisotropic Ginzburg-Landau model 

In this Appendix we provide details and further results on the derivation and the study of 
the coefficients of the anisotropic GL equations for a layered superfiuid Fermi gas, providing 
as well the equations giving the critical temperature Tc. The starting point is the Hubbard 
Hamiltonian ([2]): 

where is the fermionic operator, Uj denote the tunnelhng rates, /x is the chemical potential 
and 0" = t)i is the spin index. Furthermore i,j are site indices and N = J^ia^la^ia 
total number operator. The interaction is assumed attractive, U > 0, and we denote by n 
the total number of fermions per site {n < 2). The tunnelling rate between the sites in the 
xy planes is denoted by t\\, while the tunnelling rate along the z axis is denoted by t±. 

As discussed in the main text, in our study of the anisotropic GL and LD coefficients 
the crucial ingredient is to have the optical lattice along z and the mechanism of layering 
in the form of tunnelling between adjacent layers: with small (or absent) optical lattices in 
the x,y directions, then the anisotropic GL ([T| would be obtained. In writing (Al) we also 



that in each site of the lattice there is a total number of fermions n < 2: otherwise, one has 
a multi-band Hubbard description and again the anisotropic GL would be retrieved. 

In a path-integral setting the fermionic atoms are described by the independent Grass- 
mann fields (p and (j) [3914111 B3| HI]: the partition function is given hj Z = J T)(l)D(j)e~^f^'^''^\ 
with Sf the fermionic action in imaginary time [IQ]. First, the Hubbard-Stratonovich trans- 
formation is performed by completing the Gaussian integral: the theory has then the form 
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of a (free) fermion field interacting with the bosonic field A. The corresponding action reads 

kq k~q 



s 



u 



(A2) 



where 




kq 







(A3) 



(A4) 



lUJ 



and 



M, 



qk 



(A5) 



Afc_g 

A,_fc 

(here (3 = l/ksT, while in the main text it only denotes the coefficient of the quartic term in 



the GL free energy) . In the previous formulas eg = e?''—yU where e^""* = — 2t|| (cos fc^+cos /cj,) — 
2t_L cos kz denote the single-particle energies (in this Appendix the lattice spacing is set to 1). 
Furthermore we used the notation k = [iu; k) where the wavevectors k belong to the Brillouin 
zone [—71; vr]^ and ue = in/P are the Matsubara frequencies (with i even or odd in the case of 
bosons and fermions respectively). The sum over k corresponds to an integral over k over the 
Brillouin zone and sum over the integers i; furthermore, 0j,(r) = J2k e'^^'"-' 0,, (iw^; fc) 
{fj is the position of the site j). 
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FIG. 5. Anisotropy parameter 7^ = m±/m\i^ vs. U/D at fixed t±/t^^ = 0.5 and different particle 
densities n (from top to bottom n = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9). 
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The fermionic degrees of freedom can be integrated out obtaining the bosonic effective 
action 



The critical temperature Tc is the temperature for which the solution has Aq = where 
Ao = (A). When the interaction is weak, only the atoms around the Fermi level form 
Cooper pairs (BCS limit) [31]. In the BCS limit Cooper pairs have zero relative momentum, 
i.e. (Afe) ^ for k ^ 0, and the chemical potential is not changed from the Fermi energy 
Ep. In this limit the phase-space integration can be done by saddle point approximation 
leading to the mean-field equations. The result is that Tc ^ Tp, where Tp is the Fermi 
temperature [Ml E] . When this is not the case, one has to consider the effects of higher 
order fluctuations [IQ] and solve two coupled equations to obtain fi and Aq as a function of 
the interaction U, the hopping rates ty and t±, the particle density n and the temperature 
T. The first equation is the saddle point equation for the order parameter Aq: 




(A6) 



k 




(A7) 



where Ej^ 



e| + Aq. The other is the number equation 

n = -—F 



(A8) 
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FIG. 6. Effective anisotropy 'y'^ff a function of U/D for different anisotropics (at n = 0.5): from 
top to bottom i_LA|| = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1. 
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FIG. 7. Plot of 7^ • {t±/t\\Y vs. U/D at n = 0.5 for the same layering values t±/t\^ of Figj6| 



where (apart from a constant) 



F = ~Y. In (1 + e-%) + ^ + ^('?) 



keBZ 



(A9) 



is the free energy and is the pair polarization function, which for Aq = reads 



X 



(AlO) 



with 



To determine the critical temperature one has to solve the two Eqs. (A7)-(A8) with 



Ao = 0: in this way one obtains and /i as a function of the tunneling rates and t_L, the 



particle density n and the interaction U . We observe that the numerical solution of (A7) 



(A8) involves for the lattice case the explicit evaluation of (derivatives of) integrals over six 



variables, while in the 3D continuum case one can exploit the rotational invariance and 
actually evaluate integrals in three variables. A simplification in the numerics is obtained 
noticing that after transforming the sum over Matsubara frequencies into a contour integral 
a first order pole appears that has to be corrected by either a limiting procedure putting 
iw + O^, or by subtracting the contour integral at infinity which can be obtained analytically. 



The solution of Eqs. (A7)-(A8) for vs. U for different fillings and anisotropics is shown 



in FigjT} One also observes that, as in the continuum case, the unitary point is different 
from the maximum point. 
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FIG. 8. The different regimes spanned by a layered superfluid: (a) in the BCS regime the Cooper 
pairs are spread both along and across the planes; (b) in the Lawrence-Doniach regime the Cooper 
pairs live each in one plane - having a large in plane radius ~ and tunnel between the layers as a 
single particle; (c) in the BEC regime each pair forms a tightly bound molecule and stays mostly 
in one site. As shown in Figj4| the region where the Lawrence-Doniach model is important in the 
crossover between the BCS and BEC limits. 

If the system is near the critical temperature Tc, the average value of A is small and the 
effective action can be expanded in series [39l HS] . Stopping at the fourth order gives 

S-So = l tT{GMGM) + J2 ^ + ^ tT[{GM)% (All) 

k 

where 5*0 = lndet(G^^) — describes a non-interacting Fermi gas. The last step in finding 
a GL form for the action is to expand each coefficient around g* = (time-dependence can 
be also included by expanding in iu). The resulting action is of the same form as the GL 
free energy ([T]), with the GL coefficients m±, my, a, (3 given in Eqs. ([3])-(|4]). 

We remark that the action just derived is also the free energy when the value (A) is fixed. 
A natural question arises whether this functional corresponds to the original GL free energy, 
that is a functional of the anomalous correlation function. It can be seen that the minimum 
condition di^F[/S\ = gives {(p-^) {4>i4>^) = A/[/, so that the free energy just found gives 
the correct equilibrium values for the observables. 

Once that the GL coefficients m_|_, my, a, (3 are obtained (and their explicit evaluation 
requires the previous determination of T^. and yu), it is then possible to evaluate the anisotropy 
parameter 7 = \fm^Jrn\\. Our results for 7 as a function of the interaction f/ at a fixed 
filling n are plotted in Fig|2] We also studied the influence of the filling on the anisotropy 
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parameter: our findings are reported in FigjSj The filling modifies 7 on the BCS side, where 
the dynamics is determined by the atoms at the Fermi surface. On the other side of the 
unitary limit, the dynamics interests all the energy scales and is less susceptible to filling. 
In particular while the critical temperature for a given interaction value U increases with 
filling, the anisotropy on the BEG side remains virtually identical. 

A natural question is whether the collective behaviour enhances or not the anisotropy 
that is artificially tuned by imposing a layering through a value of smaller than one: 

for this reason we define an effective anisotropy 7e// defined through the relation 

2 2 ^-L t± / A 1 o^ 

leff = 7 ^ = T- 0^12 

tj| mj| t|| 



The reason behind the choice of the definition (A12) is that we expect that the effective mass 



along the direction a is ~ 1/ta (where a = x,y,z). Then a value of jeff larger than one 
signals an effective enhancement of the anisotropy (with respect to the "external" anisotropy 
t\\/t±). The results are drawn in Figjoj it is seen that in the BCS side the enhancement is 
more pronounced with respect to the BEC side. Again, the effective enhancement is reduced 
in the unitary limit. 

The pair tunnelling in the BEC limit is expected to approach the strong coupling value 
j_/f/ for each direction. Indeed in the BEC limit the quantity 7^ • (t±/t||)^ is found to 
approach 1 (with the convergence to the BEC limit being slower for systems with stronger 
layering), as shown in FigjTj 

Finally, in FigjS] we pictorially depict the physical regimes of a layered ultracold Fermi 
gas. 
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